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ANALYTICAL  DESCRIPTION  OF  THE  BREAKUP 
OF  LIQUID  JETS  IN  AIR 


Demetrios  T.  Papageorgiou1 
Department  of  Mathematics  and 
Center  for  Applied  Mathematics  and  Statistics 
New  Jersey  Institute  of  Technology 
Newark,  New  Jersey  07102 


ABSTRACT 

A  viscous  or  inviscid  cylindrical  jet  with  surface  tension  in  a  vacuum  tends  to  pinch  due 
to  the  mechanism  of  capillary  instability.  We  construct  similarity  solutions  which  describe 
this  phenomenon  as  a  critical  time  is  encountered,  for  two  physically  distinct  cases:  (i) 
Inviscid  jets  governed  by  the  Euler  equations,  (ii)  highly  viscous  jets  governed  by  the  Stokes 
equations.  In  both  cases  the  only  assumption  imposed  is  that  at  the  time  of  pinching  the 
jet  shape  has  a  radial  length  scale  which  is  smaller  than  the  axial  length  scale.  For  the 
inviscid  case,  we  show  that  our  solution  corresponds  exactly  to  one  member  of  the  one- 
parameter  family  of  solutions  obtained  from  slender  jet  theories  and  the  shape  of  the  jet  is 
locally  concave  at  breakup.  For  highly  viscous  jets  our  theory  predicts  local  shapes  which 
are  monotonic  increasing  or  decreasing  indicating  the  formation  of  a  mother  drop  connected 
to  the  jet  by  a  thin  fluid  tube.  This  qualitative  behavior  is  in  complete  agreement  with  both 
direct  numerical  simulations  and  experimental  observations. 
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tract  No.  NAS1-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1  Introduction 


It  is  well  known  (Rayleigh  (1878)),  that  a  circular  jet  of  finite  radius  with  a  surface  which 
supports  surface  tension  is  linearly  unstable  to  a  long-wave  capillary  instability.  Any  small 
perturbations  with  wavelengths  larger  than  the  jet  radius  grow  exponentially.  Furthermore, 
linear  theory  predicts  a  maximally  growing  wave  and  hence  a  dominant  length  scale  for 
the  instability.  According  to  linear  theory,  then,  the  dominant  wavelength  is  approximately 
9n  where  a  is  the  unperturbed  jet  radius  (this  result  is  independent  of  the  surface  tension 
coefficient).  Experiments  indicate  that  the  instability  can  lead  to  breakup,  or  pinching,  of  the 
jet  into  drops.  Clearly  the  pinching  phenomenon  is  nonlinear  since  the  initial  disturbance  has 
to  grow  to  amplitudes  of  the  order  of  the  unperturbed  jet  radius.  Linear  theory,  however, 
does  well  in  the  qualitative  prediction  of  breakup  times,  for  instance,  by  employment  of 
emphka!  arguments  such  as  e-fold  amplification  of  perturbations.  In  many  applications  the 
shape  and  jet  velocities  at  breakup  are  useful  but  cannot  be  obtained  from  linear  theory. 
This  classical  problem  has  been  studied  extensively;  experiments  have  been  carried  out  by 
Donnelly  and  Glaberson  (1966),  Cloedde  and  Yuen  (1970)  and  more  recently  Chaudhary  and 
Maxworthy  ( 1980a, b).  Weakly  nonlinear  theories  (see  below)  have  been  carried  out  by  Yuen 
(1968)  and  later  by  Chaudhary  and  Redekopp  (1980).  A  review  of  the  subject  can  be  found 
in  Bogy  (1979)  while  recent  simulations  using  boundary  integral  techniques  are  described  in 
Mansour  and  Lundgren  (1990). 

Viscosity  dominated  flows  form  a  separate  but  parallel  field.  Tomotika  (1936)  consid¬ 
ered  the  linear  stability  of  a  stationary  cylindrical  thread  of  viscous  fluid  surrounded  by  a 
second  viscous  fluid  with  surface  tension  acting  at  the  interface.  Qualitatively,  the  stability 
results  are  similar  to  inviscid  studies  with  a  maximally  growing  wave  of  the  order  of  the 
unperturbed  thread  radius.  A  more  complete  theory,  including  the  effects  of  non-uniform  jet 
velocities  can  be  found  in  Chandrasekhar  (1961),  where  it  is  shown  that  capillary  instability 
provides  linearly  growing  waves  that  scale  on  the  jet  radius.  Recently,  Tjahjadi,  Stone  and 
Ottino  (1992),  have  undertaken  an  experimental  and  numerical  study  of  breakup  of  viscous 
cylindrical  threads  of  one  fluid  in  another.  The  experiments  and  the  computations  show  that 
at  the  time  of  breakup  the  jet  tends  to  form  larger  mother  drops  joined  to  smaller  satellite 
drops  by  thin  slender  tubes.  It  is  precisely  this  regime  that  we  can  describe  theoretically 
with  excellent  qualitative  agreement  with  both  computations  and  experiments.  Besides  their 
intrinsic  interest,  such  solutions  can  be  useful  in  providing  initial  conditions  for  continuation 
of  numerical  solutions  just  before  the  change  in  topology  necessitated  by  the  pinching.  The 
analytical  description  of  local  structures  is  also  useful  in  the  determination  of  the  effect  of 
additional  physicochemical  effects  such  as  surface  active  agents  or  electrical  forces  on  the 
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pinching  process. 

The  present  approach  is  a  fully  nonlinear  one  with  interfacial  deflections  as  large  as  the 
undisturbed  jet  radius,  as  opposed  to  most  previous  weakly  nonlinear  studies.  Chaudhary 
and  Redekopp  (1980)  (see  also  Yuen  (1968))  consider  two  low  amplitude  (asymptotically 
small  but  not  infinitesimally  so)  initial  perturbations,  a  fundamental  and  a  harmonic.  These 
are  followed  up  to  cubic  order  in  the  initial  small  amplitude,  which  is  the  first  stage  when 
the  modes  interact  nonlinearly  to  produce  an  amplitude  equation.  The  methodology  is  that 
of  the  Stuart  (1960),  Watson  (1960)  classical  weakly  nonlinear  theory,  even  though  the  jet 
problem  always  has  a  band  of  unstable  waves  which  does  not  become  monochromatic  as  a 
flow  parameter  (e.g.  capillary  number  here)  is  varied.  The  results  are  therefore  expected  to 
be  valid  for  sufficiently  small  times,  but  as  Chaudhary  and  Redekopp  indicate  qualitative 
features  of  the  experiments  are  reproduced  at  times  which  seem  to  be  beyond  the  validity 
of  the  theory.  Fully  nonlinear  theories  allow  the  interfacial  amplitude  to  be  as  large  as  the 
unperturbed  jet  radius,  a  situation  which  is  essential  in  the  description  of  breakup.  This 
usually  means  that  the  problem  should  be  addressed  numerically,  and  it  is  the  objective  of 
this  work  to  present  an  analytical  alternative  to  breakup.  A  fully  nonlinear  theory  of  jet 
breakup  in  a  vacuum  was  developed  by  Ting  and  Keller  (1991),  referred  to  as  TK.  This 
is  done  by  use  of  slender  jet  theory  (see  below)  which  was  extended  by  Papageorgiou  and 
Orellana  (1993),  referred  to  as  PO,  to  describe  breakup  of  jets  of  one  fluid  into  another 
with  different  density,  with  or  without  surface  tension  and  inclusion  of  Kelvin-Helmholtz 
instability.  In  both  cases,  a  simplified  set  of  evolution  equations  involving  time  and  the  axial 
coordinate  alone,  are  derived  by  an  asymptotic  expansion  procedure  which  uses  the  ratio 
between  undisturbed  jet  radius  to  characteristic  axial  length  scale  as  a  small  parameter.  As 
noted  by  Papageorgiou  and  Orellana,  such  an  ansatz  can  be  applied  to  flows  which  have 
initial  conditions  characterized  by  a  long  axial  length  scale,  for  instance.  This  implies  that 
the  slender  jet  theory  preserves  such  geometry  throughout  the  evolution  and  at  the  time  of 
pinching  the  distance  between  successive  pinch  points  is  large  compared  to  the  undisturbed 
jet  radius.  Qualitatively,  therefore,  the  slender  jet  theory  predicts  drops  with  radii  larger 
than  the  undisturbed  jet  radius.  The  strength  of  the  slender  jet  theory,  however,  is  in  its 
prediction  of  local  structures  at  the  time  of  pinching.  As  found  in  TK  and  extended  to  other 
flows  by  PO,  a  one-parameter  family  of  similarity  solutions  is  possible.  Different  members 
of  this  family  correspond  to  local  geometrical  shapes  of  the  jet  which  are  concave,  conical 
and  convex  respectively  (see  PO  for  representative  numerical  solutions).  In  this  work  we 
show  that  if  the  slender  jet  theory  is  not  adopted  throughout  the  evolution,  the  jet  can 
pinch  according  to  a  single  member  of  the  one-parameter  slender  jet  family,  corresponding 
to  concave  shapes.  Solutions  with  zero  surface  tension  also  follow  by  simply  discarding  one 
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term  from  the  similarity  equations.  As  shown  in  PO,  the  slender  jet  equations  with  zero 
surface  tension  constitute  a  system  of  two  nonlinear  hyperbolic  conservation  laws  with  the 
final  outcome  of  the  initial  value  problem  being  the  formation  of  a  shock  after  a  finite  time. 
Clearly  the  development  of  short  axial  length  scales  is  inconsistent  with  the  fundamental 
assumptions  of  the  theory.  The  direct  theory  presented  here,  however,  does  not  assume 
the  slender  jet  approximation  throughout  the  evolution  and  so  does  not  suffer  from  the 
restrictions  outlined  above.  Consequently,  a  pinching  solution  for  zero  surface  tension  can 
be  given  for  a  flow  which  is  stable  according  to  linear  theory.  The  shapes  at  breakup  are 
monotonic  indicating  the  possibility  of  formation  of  mother  drops  joined  to  satellite  drops  by 
thin  fluid  tubes.  Such  shapes  are  seen  in  experimental  studies  of  inviscid  jets  (see  Donnelly 
and  Glaberson  (1966)). 

We  emphasize  that  slender  jet  theories  are  powerful  methods  in  the  description  of  local 
structures  in  pinching  phenomena.  In  capillary  instability  phenomena  in  three  dimensions, 
slender  jet  theories  provide  simplified  sets  of  evolution  equations  which  can  be  analyzed  for 
breakup.  Since  at  breakup  of  cylindrical  jets  the  slender  jet  assumptions  are  valid  it  is  not 
surprising  to  find  that  a  description  of  breakup  emerges.  Physically,  the  disturbance  must  be 
characterized  by  large  wavelengths  initially.  The  present  approach  corresponds  to  a  different 
physical  set-up  as  we  describe  next.  If  the  jet  is  allowed  to  undergo  natural  growth  without 
any  external  forcings  and  linear  waves  are  at  first  important,  then  linear  theory  can  be  used 
to  predict  a  maximally  growing  wave.  As  mentioned  above,  this  wave  has  wavelength  of  the 
order  of  the  undisturbed  jet  radius  and  so  when  the  evolution  enters  the  nonlinear  regime, 
the  interfacial  waves  have  axial  length  scales  comparable  to  the  undisturbed  jet  radius. 
The  slender  jet  approximation  may  not  be  appropriate,  therefore,  for  the  total  duration  of 
the  evolution  to  breakup,  even  though  it  is  appropriate  locally  at  breakup.  The  analysis 
presented  here,  assumes  that  near  the  singular  time  the  jet  radius  is  asymptotically  smaller 
than  the  axial  length  scale,  a  condition  which  may  be  inappropriate  for  the  past  history  of 
the  evolution. 

The  article  is  organized  as  follows.  Section  2  deals  with  the  breakup  of  inviscid  jets 
governed  by  the  axially  symmetric  Euler  equations.  A  system  of  differential  equations  is  de¬ 
rived  that  describes  the  local  structure  of  the  shape  and  the  potential  of  the  jet  at  breakup. 
Numerical  solutions  are  given  and  a  correspondence  is  made  with  the  slender  jet  results  of 
PO,  both  in  the  presence  and  absence  of  surface  tension.  In  Section  3  we  analyze  viscous 
axisymmetric  jets.  The  equations  studied  are  the  Stokes  equations  (creeping  motion  corre¬ 
sponding  to  small  Reynolds  numbers)  and  a  differential  equation  describing  the  jet  shape  at 
breakup  emerges.  The  analysis  leading  to  this  equation  is  more  involved  than  the  inviscid 
case  due  to  the  higher  order  of  the  field  equation  and  the  additional  stress  boundary  condi- 
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tions.  Numerical  and  analytical  aspects  of  the  equation  are  also  given.  Finally  Section  4  is 
devoted  to  discussion  and  comparison  with  experiments  and  simulations. 


2  The  theory  for  inviscid  jets 


The  equations  governing  the  fluid  motion  can  be  written  in  terms  of  a  velocity  poten¬ 
tial  4>(t,r,z)  which  is  independent  of  the  azimuthal  angle  0  for  axisymmetric  flows.  A 
system  of  cylindrical  polar  coordinates  (?-,  9 ,  z)  is  used  with  corresponding  velocity  vector 
u  =  ( u,v,w )  =  V<j).  Using  this  notation,  (f>  satisfies  Laplace’s  equation  while  at  the  free 
surface  two  conditions  are  specified,  a  kinematic  condition  and  a  normal  stresses  balance 
which  gives  the  pressure  jump  across  the  interface.  When  the  latter  condition  is  used  in  the 
2-momentum  equation  the  usual  Bernoulli  equation  arises.  Without  loss  of  generality  we 
give  the  system  for  a  static  undisturbed  jet,  noting  that  any  background  axial  velocity  can 
be  removed  by  a  Galilean  transformation.  The  equations  with  boundary  conditions  are 


<t>rr  +  ~<f>r  +  <f>zz  =  0, 


(2.1a) 


On  r  —  S(t ,  z ), 


<t>T  =  Si  +  (f>zSz, 


+  +  -  ~  (I  ”  s-  +  TTfj) (1  +  s‘‘ 


S,Y'n- 


(2.16) 

(2.1c) 


A  final  condition  which  completely  specifies  the  problem  is  regularity  of  (f>a.tr  —  0.  Equations 
(2.1a-c)  are  in  non-dimensional  form  with  the  non-dimensional  group  containing  the  surface 
tension  coefficient  (the  Weber  number)  scaled  to  unity,  in  anticipation  of  this  canonical 
rescaling  in  the  final  similarity  equations. 

In  general  (2.1a-c)  must  be  addressed  numerically.  We  construct  pinching  solutions  of 
(2.1a-c.)  by  assuming  the  following  ansatz  as  a  singular  time  ts  is  approached  from  below: 


r  =  (ts-t)ay  ,  z  =  (ts-tft  ,  S(t,z)  =  (t,-t)af(0,  (2.2  a) 


<i>{t,r,z)  =  (t,  -  tyx{i,y,  0- 


According  to  these  transformations  time  and  space  derivatives  become 

d  ,  8  |  M  9 

dt  dt  (ts  —  t)dz' 


d_ 

dr 


( ts-t ) 


d_ 

dy 


d 

dz 


{ts-t)-0 


d_ 

dy 


(2.26) 


(2.2c) 
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Physically,  we  expect  the  jet  radius  to  be  much  smaller  than  the  axial  length  scale  at  breakup. 
This  implies  that  r  <<  z  which  in  turn  yields  the  condition 

a  >  il.  (2.3) 

Inequality  (2.3)  is  the  only  necessary  assumption  for  the  theory  to  be  consistent. 

Substitution  of  the  ansatz  (2.2a-c)  into  the  Laplace  equation  (2.1a)  indicates  that  since 
a  >  /U,  the  dominant  terms  come  from  the  r- derivatives.  The  leading  order  contribution 
which  is  regular  at  the  jet  axis  y  =  0  is  independent  of  y,  and  we  must  go  to  a  higher  order 
to  obtain  a  non-trivial  balance.  This  suggests  seeking  a  solution  by  expanding  \  'a  (2.2b) 
in  appropriate  powers  of  ( t.3  —  t)  as  follows: 

4>  =  {t3-  ty  (xo(j /, o  +  {t,  -  t)0Xi(v,{)  +  {U  -  t)°X2 (y,0  +  •••)•  (2.4) 

Substitution  of  (2.4)  into  the  Laplace  equation  (2.1a)  gives 

( 1 5  -  ty~2a  ^XOvv  +  “XOyj  +  ( 1 4  -  ty+0~2a  ^Xlyy  +  “Xly^  + 

(t,  -  ty~: +  .  • .  =  0.  (2.5) 

It  can  be  seen  from  (2.5)  that  the  first  stage  where  a  regular  y-dependent  solution  is  possible, 
arises  if  the  Xi  terms  are  in  balance  with  the  Xo«  term.  This  yields 

a  =  p.  (2.6) 

The  solutions  for  Xo  a°d  Xi  which  are  regular  at  y  =  0  are 

Xo  =  4(0,  XI  =  -  +  B{  (),  (2.7a,  b) 

where  A  and  B  are  unknown  functions  of 

Next  we  consider  the  kinematic  condition  (2.1b).  Substitution  of  the  ansatz  (2.2)  and 
(2.4)  along  with  the  solutions  (2.7a,b)  yields  the  equation 

<(,  -  <r+s-“  (- i/4«)  =  — Q(( .  -  ()“-'/ 

+D(t,  -  +  ((,  -  t),+“-M4{/ (  +  . . .  .  (2.8) 

A  few  comments  about  (2.8)  are  in  order.  First,  the  two  terms  which  contain  A  are  multiplied 
by  the  powers  7  +  /?  —  a  and  7  4-  a  —  20  of  (t3  —  t )  respectively.  Since  a  —  |/3  from  (2.6)  these 
terms  are  in  balance  and  provide  no  new  information.  An  expression  for  7  is  next  found  by 
making  a  balance  of  all  four  leading  order  terms  in  (2.8).  This  gives 

7  =  2/3-1, 


(2.9) 


and  the  kinematic  condition  becomes  to  leading  order 

+  W  +  (2.10) 

where  primes  denote  ^-derivatives. 

Another  equation,  as  well  as  determination  of  the  free  parameter  [i  is  found  from  the 
Bernoulli  equation  (2  lc).  Substitution  of  the  ansatz  (2.2),  (2.4)  and  retention  of  the  largest 
possible  terms  as  t.  — >  ts— ,  gives  the  following  equation: 

-7 (t,  -  ty-'A  +  (t,  -  +  (ts  -  =  ~(ts  -  tyaf-\  (2.ii) 

We  arrive  at  (2.11)  above  by  dropping  any  terms  that  are  smaller  than  the  ones  retained 
using  the  fact  that  a  >  ft  alone.  Next,  we  balance  the  different  terms  in  the  equation.  The 
terms  of  orders  ( t ,  —  t)y~l  and  ( t3  —  t)2'y_2/3  are  in  balance  by  virtue  of  (2.9)  above.  The 
surface  tension  term  containing  /-1  is  chosen  to  balance  the  other  terms  and  produces  a 
coupling  with  (2.10)  between  inertial  and  capillary  forces.  This  is  possible  if 

7  —  1  =  —or. 


Use  of  (2.6)  and  (2.9)  that  provide  a  and  7  in  terms  of  /i  yields  a  value  for  (3  and  hence  or 
and  7  follow, 


6  4  1 
<*=-  ,  p  =  —  ,  7  =  — , 

(2.12) 

while  equation  (2.11)  becomes 

1  .  4,  .  1  1 

— -4+-{  At  +  -A(--j. 

(2.13) 

Equations  (2.10)  and  (2.13)  form  a  coupled  system  which  must  in  general  be  addressed 
numerically.  Before  presentation  of  such  results  we  cast  (2.10),  (2.13)  into  a  form  which  is 
readily  comparable  with  the  one-parameter  family  equations  of  PO.  We  define  G  =  A ^  and 
/  =  \/F,  differentiate  (2.13)  with  respect  £  to  yield  the  system 

-jF+4-(F'  +  (GF)'  =  0,  (2.14o) 

\g  +  *-i(f  +  \(G‘)'  =  -(F-,/2)'.  (2.146) 

Equations  (2.14a,b)  are  exactly  the  equations  found  from  slender  jet  theory  corresponding 
to  ()  =  |.  In  the  absence  of  surface  tension  the  term  on  the  right  hand  side  of  (2.14b)  is 
dropped.  As  shown  below,  zero  surface  tension  solutions  are  possible  even  though  the  flow  is 
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linearly  stable.  Before  presentation  of  numerical  solutions  we  give  the  asymptotic  behavior 
of  (2.14a,b)  for  large  |£|  (for  more  details  see  PO).  In  the  absence  of  surface  tension  we  find 


F~£3  G~C*  £  —  oo,  (2.15a) 

F~(“0-1  G  =  £  — ►  — oo.  (2.156) 

Note  that  the  shapes  predicted  for  zero  surface  tension  grow  on  one  side  (as  £  — ♦  oo  here) 
and  decay  at  the  other  end  of  the  domain.  This  is  in  contrast  to  the  surface  tension  case  as 
the  following  asymptotic  forms  indicate, 

F  ~  £3  G~C*  £-»oo,  (2.16a) 

F~(-03  G~  (-O'*  ( -  -oo.  (2.166) 

Numerical  solutions  of  (2.14a,b)  showing  the  variation  of  F  and  G  with  £  are  given  in 
Figures  1  and  2  for  zero  surface  tension.  Verification  of  the  asymptotic  forms  (2.15a,b)  is 
included  in  Figures  3  and  4  for  the  behavior  of  F  and  G  respectively.  This  is  achieved  by 
constructing  logarithmic  plots  of  the  numerical  solutions  for  large  positive  and  negative  <f. 
The  slopes  of  the  lines  so  obtained  were  computed  from  the  data  by  a  least  squares  fit  and  are 
indicated  on  the  Figures.  It  is  seen  that  agreement  is  very  good  and  improves  as  the  range  of 
integration  is  extended.  The  results  for  non-zero  surface  tension  are  presented  in  Figures  5 
and  6  with  the  corresponding  logarithmic  plots  confirming  the  asymptotic  behavior  in  Figures 
7  and  8.  The  main  qualitative  difference  between  the  two  cases  is  that  in  the  absence  of 
surface  tension  the  shape  is  monotonic  tending  to  zero  at  one  end  of  the  similarity  region 
and  growing  algebraically  at  the  other  end,  while  in  the  presence  of  surface  tension  the  shape 
function  grows  algebraically  (at  the  same  asymptotic  rate)  at  both  ends  of  the  domain.  The 
former  type  of  behavior  is  also  obtained  for  highly  viscous  jets  analyzed  in  Section  3. 


2.1  Conservation  integrals  at  breakup. 


The  inviscid  jet  flow  is  conservative  and  solutions  to  the  problem  are  required  to  conserve 
the  total  momentum  and  energy  as  well  as  mass.  It  can  be  seen  from  the  solutions  just 
constructed  that  velocities  are  becoming  unbounded  as  the  jet  radius  shrinks  to  zero.  It  is 
important,  therefore,  to  verify  the  consistency  of  these  solutions  with  conserved  quantities  of 
the  system.  We  define  the  total  mass,  momentum  and  energy  of  an  axially  periodic  portion 
of  the  jet  by  /,,  /2  and  I3  respectively.  These  are  given  by 


h 


(2.17) 
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(2. IS) 


I2  =  2np  f  f  (<?i  +  0i)l/'*rdrd:, 

J-P  J u 

rl>  rS  J  ft) 

l\  =  27rp  /  /  -'Vr  +  $i)rdvdz  +  27rc  /  .S'dr, 

J-[)J  o  * 


(2.  IS) 


where  <r  is  the  surface  tension  ro  ,4  lent  and  the  quantities  (2. 1 7)-(2. 19)  are  in  dimensional 
forms.  I  he  contributions  of  the  integrals  /,.  and  /: j  in  the  similarity  region  are  found 
by  substitution  of  the  h-ading  order  solutions  from  the  ansatz  (2.2a-c)  to  give  the  following 
quantities 


l\  ~  ~  0  7  /  /•’</£, 

J  —  OC* 

/,  ~  (<3  -  /)t  y  x  F(;d{. 
/« ~  (/,  -  /)*  r  f<7^. 


(2.19) 

(2.20) 

(2.21) 


The  contributions  to  the  energy  integral  (2. IS)  of  the  kinetic  and  surface  energy  in  the 
similarity  region  are  of  the  same  order  and  form  and  both  appear  in  (2.21);  the  expression, 
therefore,  is  valid  both  in  the  absence  and  presence  of  surface  tension.  We  have  shown 
analytically  and  numerically  that  F  and  (1  become  unbounded  as  |£|  — >  cc.  In  order  to 
evaluate  the  contributions  (2. 1 9)- (2.2 1 )  on  the  large  scale  flow  we  need  to  transform  to  outer 
variables.  The  main  contributions  to  the  integrals  come  from  the  asymptotic  forms  (2.15a,b) 
for  zero  surface  tension  or  (2.16a,b)  for  non-zero  tension.  The  behavior  at  infinity,  however, 
is  the  same  for  both  cases  and  the  results  that  follow  are  equally  valid.  Consideration  of  the 
leading  order  contributions  of  (2. 19)-(2.21 )  gives 


/i  ~  (*,  -  t)T£*  =  (2.22) 

/,  ~  (t,  -  t)T£T  =  (t,  -  t)%z'r,  (2.23) 

/.)  ~  (t,  ~  (2.24) 

We  can  conclude  from  (2.22)-(2.24)  that  as  the  singular  time  is  approached,  the  solutions 
constructed  here  are  consistent  with  the  conservation  integrals  of  the  problem.  More  specif¬ 
ically  the  total  mass,  momentum  and  energy  in  the  breakup  region  tends  to  zero  as  the 
singular  time  is  reached.  This  theoretical  result  can  be  of  used  as  a  criterion  in  numeri¬ 
cal  simulations,  for  example,  since  near  the  singular  time  the  total  energy  of  the  system  is 
concentrated  away  from  the  similarity  region. 


3  Theory  for  viscous  jets 

In  this  section  we  consider  the  collapse  of  a  viscous  thread  of  fluid  of  undisturbed  radius  H 
under  capillary  instability.  In  numerous  applications  characteristic  Reynolds  numbers  are 


cS 


small  and  the  flow  is  governed  by  the  Stokes  equations.  Such  flows  are  driven  by  capillary 
forces  which  in  turn  provide  a  scale  for  the  flow  velocities.  In  non-dimensionalizing  the 
equations,  therefore,  the  following  scales  are  used: 

a  a  uR _ 

(*',-)  —  R(r,  c )  (u,w)=-{u,w)  p=- -p  t=  - /. 

p  R  a 

where  a  is  the  surface  tension  coefficient  and  p  the  fluid  viscosity.  The  non-dimensional 
equations  and  interfacial  boundary  conditions  become  (dropping  the  bars): 

A  u - -u  =  pr,  (3.1  a) 

r* 

A  w  =  p:,  (3.16) 

-  {vu)r  +  w:  =  0,  (3.1c) 

r 

A  -  —  — 

dr2  v  dr  +  dz2 

On  r  —  S(t,  ;)  we  have 

(u:  +  ay )( 1  -  S'i)  +  2 urSz  -  ‘2wzS:  =  0,  (3. Id) 

p  —  2 ur  -  (—  p  +  2w.)S'*  +  2(u-  +  n\)$.  = 

+.^))(1+  5^)-,/2,  (3.1c) 

u  =  $t  +  w$z.  (3.1/) 

The  interfacial  conditions  (3. 1  d-f )  represent  the  tangential  stress  balance,  normal  stress 
balance  and  the  kinematic  condition  respectively.  An  additional  condition  is  regularity  of 
the  velocity  field  on  the  jet  axis  r  =  0.  The  problem  ( 3 . 1  a- f )  poses  a  formidable  analytical 
task.  Analytical  studies  are  usually  confined  to  linear  stability  and  in  the  fully  nonlinear 
regime  the  problem  has  been  addressed  numerically  (see  Introduction).  In  order  to  draw' 
up  an  analogy  with  the  inviscid  analysis  of  Section  2,  we  can  introduce  a  streamfunction  i!> 
defined  by  u  =  —  J-i />-,  w  =  ^r/y  so  that  the  continuity  equation  (3.1c)  is  satisfied.  Elimination 
of  p  between  (3.1a,b)  yields  a  single  equation  for  ?/y  namely 

EAx[>  =  0.  E2  =  A--—.  (3.2) 

r  dr 

The  flow  field  in  the  viscous  case  is  governed  by  a  fourth  order  equation  as  opposed  to  the 
second  order  Laplace  equation  for  the  potential  in  inviscid  flows.  As  seen  below  the  analysis 
of  viscous  jet  pinching  is  slightly  more  involved.  In  what  follows  we  choose  to  work  with 
primitive  variables  rather  than  with  the  streamfunction  ?/.>. 
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Following  the  ideas  developed  in  Section  2,  we  look  for  singular  terminal  states  of  (3.1a-f) 
according  to  the  ansatz 

r  =  (t,  ~  try  Z  =  (ts  -  t)^  S  =  (t,  -  0“/(0, 

W  =  (t,  -  tyw(t,  y,  0  u  =  (ts  -  r+a-0lJ(t.,  y,  0  v  =  (ts  -  t)-aP.  (3.3) 

The  expression  for  u  follows  from  the  continuity  equation  (3.1c)  once  w  is  specified,  and  the 
scaling  for  the  pressure  p  is  a  consequence  of  the  normal  stress  balance  equation  (3.1e)  since 
capillary  instability  drives  the  dynamics.  As  before,  the  geometry  of  pinching  dictates  the 
inequality 

a  >  fi.  (3.4) 

Consideration  of  equation  (3.1b)  for  w  indicates  that  if  a  leading  order  balance  is  made 
between  the  radial  derivative  terms  of  w  and  the  pressure  gradient  term  (this  requires  7  — 
o  —  fi),  the  following  lending  order  solution  for  W  arises 

W  =  \v2P(  +  ^«),  (3.5) 

where  A(£)  is  an  unknown  function.  An  inconsistency  appears  now  if  (3.5)  is  substituted 
into  the  tangential  stress  balance  equation  (3. Id).  The  leading  order  contribution  to  (3. Id) 
after  substitution  of  the  scalings  (3.3)  is  simply  Wy  =  0  on  y  =  /(£)•  Application  of  this 
condition  to  the  solution  (3.5)  found  above  implies  that  P^f(i)  =  0  leading  to  an  inconsistent 
i  ution  P  =  const..  A  consistent  solution  is  constructed  by  expanding  W  (and  therefore  U) 
in  appropriate  powers  of  (t3  —  t )  with  the  leading  order  contribution  of  W  being  independent 
of  y  so  that  the  tangential  stress  balance  is  satisfied  identically  to  leading  order.  Higher  order 
corrections  from  (3. Id)  then  enter  to  force  a  nontrivial  solution.  The  appropriate  expansions 
for  w  and  u  become,  then, 

w  =  (ts-  ty  (w0(0  +  (t,  -  tfw,  +  (ts  -  tyw2  +  ...),  (3.6a) 

u  =  (ts-  ty+a-e  (u0(y,  o  +  (t,  -  tfu ,  +  (tt  -  tyu2  +  ...).  (3.66) 

The  solution  for  Uo  which  is  regular  at  y  =  0  follows  from  (3.1c) 

U0  =  ~\yWoi.  (3.7) 

Next  we  obtain  the  solutions  for  W\  and  U\.  Using  (3.3)  and  (3.6a)  in  equation  (3.1b)  gives 

('•  -  <r-!"  (|j  + l-^)  («.  -  +  (<.  -  trw,  +  ...) 

+(ts  -  ty~wWoa  +  ...  =  (ts-  t)-a-VPt.  (3.8) 
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Balance  of  the  leading  order  terms  involving  W0,  W\  and  P  gives 

a  =  \&  1  =  (3.9) 

and  a  solution  for  W\  follows  by  integration  with  respect  to  y  since  the  pressure  is  independent 
of  y  by  the  normal  stress  balance  (see  comments  below  also) 

w,  =  l-y\P(  -  WHi)  +  A((),  (3.10) 

where  A  is  some  function  of  (.  Using  (3.10)  and  (3.1c.)  gives 

Vi  =  -  Wotu)  +  \yAt.  (3.11) 

Implicit  in  the  solutions  (3.10),  (3.11)  is  the  assumption  that  Py  =  0.  Using  the  values  (3.9) 
for  a  and  7  shows  that  the  leading  order  balance  in  (3.1a)  is 


\d_ 
dy2  +  y  dy 


(3.12) 


Now  substitution  of  the  solution  (3.7)  for  U0  into  (3.12)  shows  that  the  left  hand  side  is  zero 
and  so 

Py  =  0.  (3.13) 

Equation  (3.13)  is  crucial  for  the  consistency  of  the  theory. 

Next  we  consider  the  tangential  stress  balance.  Using  the  values  (3.9)  gives  the  following 
equation  to  leading  order 


Vot  +  Wiy  +  2 Uoyf  -  2 Wotf  =  0  y  =  f(().  (3.14) 

Use  of  solutions  (3.7)  and  (3.10)  into  (3.14)  gives  an  ordinary  differential  equation  for  Wo(£)i 

{fK)'  =  \fr'-  (3.15) 

An  expression  for  P  in  terms  of  W0  is  available  from  the  normal  stress  balance  equation.  It 
is  found  that  to  leading  order  several  terms  in  this  equation  are  in  balance  (according  to  the 
previously  found  values  (3.9))  the  result  being 

p  -  2Uoy  =  J  y  =  m- 

Substitution  of  this  into  (3.15)  and  integration  gives  the  solution 

Wi  =  -^J  +  J2,  (3.16) 
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where  k  is  a  constant  of  integration.  With  Wo  known  in  terms  of  /,  an  equation  for  /  arises 
from  the  kinematic  condition  (3. If)  and  use  of  the  ansatz  (3.3)  along  with  the  appropriate 
derivative  transformations  (see  Section  2  also).  To  leading  order,  therefore,  the  kinematic 
condition  becomes 

(ts  -  tr+a-fiUo  =  (/.  -  +  HU')  +  (ts  -  t)MW0f. 

A  balance  of  terms  is  achieved  if  7  +  0  —  /3  =  c*  —  1  which  along  with  (3.9)  gives 

«=1  H  =  l  7  =  “-  (3.17) 

The  kinematic  equation  is 

(Wo+|{)/+(-l  +  iwj)/  =  0, 
which  on  elimination  of  Wq  from  (3.16)  becomes 

(*.  + ft)  r+(-/- * +  £)-«•  (»■* ») 

Equations  (3.16)  and  (3.18),  therefore,  provide  a  coupled  system  of  first  order  equations  to 
be  solved  subject  to  conditions  at  £  =  0,  for  instance.  In  subsections  3.1  and  3.2  that  follow 
we  consider  analytical  and  computational  aspects  of  this  system. 

First,  some  comments  on  the  consistency  of  our  theory  with  the  original  assumption  that 
the  flow  is  in  the  Stokes  regime  are  necessary.  It  is  seen  from  the  values  (3.17)  that  the 
axial  and  radial  velocities  (see  (3.6a,b))  have  size  0((ts  —  t)~1^3)  and  0(1)  respectively  at 
the  time  of  pinching.  The  magnitude  of  the  velocity  vector  is  becoming  infinite,  therefore, 
and  the  possibility  of  violation  of  the  Stokes  flow  assumption  arises.  The  terms  ignored  are 
the  unsteady  and  convective  terms  of  the  Navier-Stokes  equations  which  read 

u 

Re(ut  +  UUT  +  UWZ)  =  ~pT  +  A  U - (NS\) 

rz 

Re(wt  +  UWr  +  WWZ)  —  —pz  +  A  W,  (N  52) 

where  is  the  Reynolds  number.  In  deriving  the  similarity  equations  (3.16),  (3.18)  we 

use  information  from  Uo,U\  in  (NSl)  and  Wo,W\  in  (NS2).  The  smallest  terms  considered, 
therefore,  are  of  0((t,s  —  t.)~4/ 3)  in  (NSl)  and  0((<s  —  f)-5/3)  in  (NS2),  coming  from  the 
Laplacian  of  if]  and  W]  respectively.  The  largest  possible  terms  which  would  enter  through 
the  unsteady  and  convective  terms  on  the  left  hand  side  of  (NSl)  and  (NS2)  are  easily 
calculated  to  be  of  orders  0(((s  —  t)~x)  and  0((ts  —  t)~4^3)  respectively.  As  t  — >  ts— , 
therefore,  the  validity  of  the  Stokes  approximation  and  the  ansatz  leading  to  (3.16),  (3.18) 
is  consistent;  higher  order  terms  eventually  enter  via  unsteadiness  and  nonlinearity  in  (NSl ) 
and  (NS2)  but  are  of  no  consequence  to  the  leading  order  balances  constructed  here. 
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3.1  Analysis  of  the  similarity  equations 

Elimination  of  W0  between  (3.16)  and  (3.18)  yields  a  single  second  order  equation  for  / 
which  can  be  written  as: 


1 


/+-/-- 


k\d2f  fP  +  f-f 


2  de 


3/ 


'df 

A, 


(3.19) 


At  first  sight,  equation  (3.19)  appears  to  admit  a  class  of  even  solutions  which,  without  loss 
of  generality,  may  satisfy  /( 0)  =  /o,  /'( 0)  =  0.  A  series  solution  for  small  4,  however,  with 
the  above  initial  conditions  is  found  to  yield  a  constant  solution  /(£)  =  /0.  For  general 
values  of  4  it  can  be  shown  that  acceptable  solutions  to  (3.19)  produce  monotonic  functions. 
Two  distinct  cases  need  to  be  considered,  (a)  k  <  0,  (b)  k  >  0. 

We  consider  the  case  k  <  0  first.  If  we  assume  that  f  vanishes  at  some  finite  value 
of  4  —  4o  say,  then  by  (3.18)  we  deduce  that  Wq(4o)  =  2.  Using  this  in  (3.16),  however, 
requires  +  p(t0)  =  ^  which  >s  impossible  since  k  <  0  and  /(4o)  >  0,  showing  that  /' 

is  non- vanishing  and  so  /(4)  is  monotonic.  This  result  can  also  be  proved  by  construction 
of  closed  form  expressions  for  /'.  We  integrate  (3.19)  once  with  respect  to  4  so  that 

df 


d( 


=  AG(f), 


d  n  <nt  r\\\  -  f2  +  f~  2k 

a/  f  3f(p  +  if-±y 


(3.20) 


6J  2  > 

where  A  is  a  constant.  Equation  (3.20)  can  be  integrated  explicitly.  There  are  three  distinct 
cases:  (i)  k  —  0,  (ii)  k  <  —  T  (iii)  —  T  <  k  <  0,  giving  the  following  solutions. 


Case  (i),  k  =  0 


G(f)  = 


P 


(/+ l/6)t 


(3.21a) 


Case  (ii).  k  <  -1/72  (k  =  -2A2) 


<?(/)  = 


/ 


exp 


(/*  +  }/  + A2)*  W*-* 


Tan 


(3.216) 


Case  (iii),  -1/72  <  k  <  0  (//  =  (£  +  t^)1/2) 


3  "t  9#j 


gw = /</ + ^ + + ^  -  A 


(3.21c) 


These  expressions  show  that  G(f),  and  hence  /',  is  non-vanishing;  cases  (i)  and  (ii)  are 
straightforward  while  in  case  (iii)  the  observation  |//|  <  yields  the  desired  conclusion. 
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When  k  >  0  the  expression  for  G  is  identical  to  (3.21c)  with  \i  as  given  there.  It  follows 
that  for  k  >  0, 

H  <  "  <  °°’ 


with  equality  on  the  lower  bound  already  covered  by  the  case  k  =  0  in  (3.21a).  Using  this 
inequality  we  see  from  (3.21c)  that  at  positions  where  /  =  /x  —  the  derivative  /'  can  (i) 
become  infinite,  (ii)  be  finite  and  non-zero,  (iii)  become  zero,  depending  on  the  value  of  the 
exponent  —  |  In  terms  of  /x  these  three  cases  are 


1  1  1 
tl>t  ’  '‘=3  ’  >‘<r 

respectively.  In  what  follows  we  show  (by  sketching  the  behavior  rather  than  presenting 
rigorous  proofs)  that  case  (ii)  alone  allows  admissible  solutions. 

We  consider  the  range  |  <  /x  <  oo  first.  As  mentioned  above  the  slope  /'  becomes 
infinite  at  positions  £o  given  implicitly  by  /(£ 0)  =  /x  —  The  interfacial  slope  becomes 
double  valued;  a  local  analysis  in  the  vicinity  of  /o  =  /x  —  ^  gives 

X  -  6)  -  (/ - /«)’  J  <  9  <  2. 

We  note  that  an  infinite  slope  in  /  violates  the  ansatz  established  by  the  scalings  (3.17)  since 
spatial  variations  in  the  axial  direction  become  much  larger  than  the  variation  in  the  radial 
direction. 

Next,  we  consider  the  values  ^  <  /x  <  |.  It  is  seen  from  (3.21c)  that  the  derivative 
vanishes  when  /  =  ^  ^  since  the  exponent  of  the  last  factor  satisfies  0<— |-fi^<l. 

This  inequality  for  the  exponent,  however,  implies  that  higher  derivatives  of  the  interfacial 
amplitude  in  the  similarity  region  are  unbounded  and  this  case  is  dropped  since  we  are 
interested  in  smooth  solutions. 

The  final  possibility  is  \i  =  |  (k  =  ^)  which  casts  (3.21c)  into 

G(f)=f(f  +  ^ri-  (3.21(f) 

Equations  (3.21a)  and  (3.2 Id)  can  be  integrated  to  yield  closed  form  expressions  for  £  in 
terms  of  /;  the  solutions  for  /  are  implicit,  therefore,  and  are  not  very  illuminating  besides 
the  checks  they  provide  for  the  numerical  and  asymptotic  work.  These  solutions  are 
k  =  0,  (3.21a) 


k  =  1,  (3.21  d) 


if+ll 


f 


+  3/(/;  6i/3 


)  —  ^  B0, 


(3.22  a) 


14 


(3.226) 


where 


51/3 

1(f-,T¥r3)  =  M  +  Bu 


/(/;  ic)  =  ^(/  +  «3)*  +  *2/»  ((/  +  /c3)’  -  k) 


—n2ln 


((/  +  *3)> 


+  «(/  T  K3)^  +  K2 


+ 


4/c 


sfi 


tan 


-l 


2  (/  +  /c3)*  +  /c 
/c\/3 


The  constants  A0,Ai,Bo,fii  are  found  from  the  initial  conditions  on  /  and  /'  at  (  =  0  for 
instance.  The  behavior  of  the  solutions  (3.22a,b)  for  large  |£|  are  easily  obtainable.  Without 
loss  of  generality  we  assume  Ao,i  >  0  to  obtain 

lb  =  0 


/~(^)V 

i  -» 00, 

f — Li  ^ 

— »  —00. 

(3.23a,  6) 

/~(^)3v 

£  -*•  00, 

/~exP(  51/3  0 

£  — ►  —00. 

(3.24a,  6) 

For  general  values  of  k  the  asymptotic  behavior  is  more  easily  obtainable  directly  from 
equation  (3.19).  Assuming  that  /'( 0)  >  0,  we  look  for  solutions  that  grow  as  £  — ►  +oo  but 
decay  as  £  — >  —  oo  (such  behavior  is  supported  by  the  closed  form  solutions  (3.22a,b)  as  well 
as  numerical  calculations).  To  leading  order,  equation  (3.19)  takes  the  forms 

J  =  \-j  /»'•  (325a) 

j  =  j(k?0)  ,  ~  ~  2jr  (k  —  0)  /«  1.  (3.25i.,c) 

It  can  be  seen  that  for  large  positive  £  the  asymptotic  form  is  independent  of  k,  in  agreement 
with  (3.23a)  (3.24a)  above.  The  asymptotic  forms  are,  therefore, 

/  ~  £3/2  £  -+  oo,  (3.26a) 

/  ~  exp(— q|£|)  (fc^O)  (*  =  0)  (  -»  -oo,  (3.266) 

where  q  is  a  positive  constant  which  can  be  determined  numerically.  Given  the  behavior  of 
/  the  asymptotic  forms  for  the  axial  jet  velocity  Wo  follow  from  equation  (3.16), 

W0  ~  const.  +  0{CV2)  £  - ►  oo,  (3.27a) 
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W0  ~  exp(2q\(\)  (k  jk  0)  W0  ~  |£|3  (*  =  0)  -oo.  (3.276) 

In  the  following  section  we  present  numerical  solutions  of  the  similarity  equations  with 
particular  emphasis  on  the  verification  of  the  analytical  findings  of  this  section. 

3.2  Numerical  solutions  of  the  similarity  equations 

Equations  (3.16),  (3.18)  or  equivalently  equation  (3.19)  were  solved  numerically  by  specifying 
conditions  at  £  =  0  and  integrating  out  to  +oo  and  — oo.  Typically  we  worked  with  the  second 
order  equation  (3.19)  with  initial  conditions  /( 0)  =  /o  and  /'( 0)  =  /i  >  0.  Integration  to  an 
appropriately  large  positive  and  negative  value  of  £  was  achieved  by  a  fourth  order  Runge- 
Kutta  method.  Representative  solutions  are  considered  for  the  three  cases  (i)  k  <  0,  (ii) 
k  —  0,  (iii)  k  =  We  note  that  the  distinctions  between  different  negative  values  of  k  are 
not  important  in  the  numerical  calculations. 

Figure  9  shows  the  results  for  /  and  Wo  of  a  case  having  k  =  —  1  and  initial  conditions 
/( 0)  =  1  and  /'( 0)  =  1.  To  establish  the  asymptotic  behavior  of  the  solutions  for  large 
|£|  we  construct  logarithmic  plots  of  the  relevant  parts  of  the  solutions.  Results  for  large 
negative  £  are  shown  in  Figure  10  which  shows  the  variation  of  ln(f)  and  ln(w)  with  The 
computed  slopes  included  on  the  curves  were  calculated  by  a  least  squares  fitting  of  the  data. 
According  to  the  asymptotic  forms  (3.26b)  and  (3.27b),  w  grows  exponentially  at  twice  the 
rate  at  which  /  decays  exponentially.  The  computed  slopes  of  2.269  and  —1.138  for  ln(w) 
and  ln(f)  respectively  are  in  full  agreement  with  this  finding.  The  asymptotic  form  (3.26a) 
for  large  positive  £  is  confirmed  in  Figure  11  with  very  good  agreement.  Verification  of  the 
behavior  of  W0  was  also  obtained  by  plotting  the  variation  of  W0  with  4-1'2  for  large  £. 

For  A:  =  0  and  the  same  initial  conditions  as  above  qualitatively  similar  solutions  are 
obtained  for  both  /  and  W0  and  representative  solutions  are  given  in  Figure  12.  A  comparison 
between  the  numerical  and  asymptotic  solutions  for  large  positive  £  according  to  (3.23a)  is 
described  next.  The  behavior  (3.23a)  is  equivalent  to 

HJ)  ~  |h  (|)  +  I'n  Q  + 

The  value  of  the  intercept  follows  from  (3.20)  when  the  values  /( 0)  =  /'( 0)  =  1  are  used, 
and  is  equal  to  —0.224.  A  logarithmic  plot  of  the  numerical  solution  produced  a  slope  of 
1.498  and  an  intercept  —0.208  (these  results  were  obtained  by  integration  to  £  =  500;  these 
values  improve  as  the  maximum  integration  range  is  increased).  The  results  for  k  =  5  are 
qualitatively  similar  to  those  described  above  and  are  not  included  here. 
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4  Conclusions 

A  formal  theory  has  been  developed  to  describe  the  breakup  process  of  highly  viscous  or 
inviscid  jets  in  air.  The  similarity  solutions  obtained  here  constitute  inner  solutions  in  the 
vicinity  of  regions  where  the  jet  radius  is  going  to  zero  which  when  matched  with  outer  solu¬ 
tions  of  the  Stokes  or  Euler  equations  respectively,  provide  a  full  description  of  the  breakup 
phenomenon.  In  general  the  outer  solutions  are  fully  nonlinear  and  must  be  computed  numer¬ 
ically.  In  numerical  calculations  of  jet  breakup  (by  boundary  integral  methods  for  example), 
it  is  important  to  have  a  rational  way  of  continuing  the  computations  beyond  the  change  in 
topology  necessitated  by  the  physics.  The  similarity  solutions  presented  here  provide  such 
a  possibility  since  accurate  initial  conditions  for  the  pinched  part  of  the  evolution,  consis¬ 
tent  with  the  full  equations  of  motion  can  be  easily  calculated.  Direct  simulations  near  the 
pinching  times  contain  stiffness  problems  which  are  likely  to  produce  inaccuracies.  Another 
example  of  the  need  of  accurate  initial  conditions  for  the  jet  evolution  beyond  breakup  is  the 
experimentally  and  computationally  observed  behavior  of  a  residual  fluid  tube  at  breakup 
being  accelerated  by  surface  tension  to  form  a  jet  which  can  go  through  the  mother  drop 
it  was  connected  to.  The  speeds  of  such  jets  depend  on  interfacial  tension  forces  which  are 
proportional  to  the  interfacial  surface  area  and  hence  the  shape  at  breakup. 

As  an  illustration  of  the  mechanics  of  the  matching  process  we  consider  the  results  of 
Section  3  for  highly  viscous  jets.  Matching  is  achieved  as  |£|  — ►  oo  and  so  the  asymptotic 
forms  of  the  solutions  are  essential  in  achieving  this.  Using  (3.26a),  for  example,  we  have 
/  ~  £3/2  as  £  — ►  oo.  Changing  to  physical  (outer)  variables  using  the  ansatz  (3.3)  and  the 
scales  (3.17)  gives  the  following  condition  to  be  satisfied  by  the  outer  solution  as  z  — ►  0  (we 
have  assumed,  without  loss  of  generality  that  pinching  occurs  at  z  =  0), 

S  ~  Z2. 

The  multiplicative  constant  that  goes  along  with  this  behavior  and  which  can  be  found 
from  the  outer  solution,  fixes  the  inner  problem  whose  solution  is  easily  obtained  from  the 
similarity  equations. 

We  have  also  found  a  qualitative  difference  between  the  breakup  forms  of  inviscid  jets 
with  surface  tension  and  their  highly  viscous  counterparts.  In  the  former  case  our  solutions 
indicate  that  at  breakup  the  interface  is  more  or  less  symmetric  with  growth  to  mother  drops 
on  either  side  of  the  pinch  point.  In  the  latter  flows  (as  well  as  inviscid  jets  with  zero  surface 
tension),  our  similarity  solutions  indicate  that  at  breakup  there  are  thin,  and  relatively  long, 
fluid  tubes  which  taper  to  zero  diameter  on  one  side  and  grow  to  connect  to  a  mother  drop 
on  the  other  side.  These  findings  are  in  full  qualitative  agreement  with  photographs  from 
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various  experiments  on  both  highly  viscous  jets  (Tjahjadi  et  al.  1992)  as  well  as  inviscid 
ones  (Donnelly  and  Glaberson  1966,  Chaudhary  and  Maxworthy  1980a, b). 
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•  Figure  1  Inviscid  jet  with  zero  surface  tension;  similarity  solution  /(£)• 
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•  Figure  2  Inviscid  jet  with  zero  surface  tension;  similarity  solution  g(£). 
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Asymptotic  benavior  of  t  at  -.-infinity 


•  Figure  3  Inviscid  jet  with  zero  surface  tension;  asymptotic  behavior  of  /(£). 


Asymptotic  bahavior  of  g  at  infinity 


•  Figure  4  Inviscid  jet  with  zero  surface  tension;  asymptotic  behavior  of  g({). 
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•  Figure  5  Inviscid  jet  with  surface  tension;  similarity  solution  /(£) 


Velocity  in  self-similar  region 
beta=4 n 
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•  Figure  6  Inviscid  jet  with  surface  tension;  similarity  solution  </(£). 
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Behavior  at  plus  infinity 
k-i  f(o>-i  r(o>-i 


stape»i.49 


Figure  11  Highly  viscous  jet  with  surface  tension.  Asymptotic  behavior  as  £ 
k  =  -1,  initial  conditions  /(0)  =  1,  /'( 0)  =  1. 


Similarity  solutions 
k-0  f(0)-1  f(0)-1 


Figure  12  Highly  viscous  jet  with  surface  tension.  Similarity  solutions  W0(0  /(£). 
k  —  0,  initial  conditions  /( 0)  =  1,  /'( 0)  =  1. 
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